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Abstract. A continuous Markovian model for truncated Levy random walks is proposed. It generalizes 
the approach developed previously by Lubashevsky et al. Phys. Rev. E 79, 011110 (2009); 80, 031148 
(2009), Eur. Phys. J. B 78, 207 (2010) allowing for nonlinear friction in wondering particle motion and 
saturation of the noise intensity depending on the particle velocity. Both the effects have own reason to 
be considered and individually give rise to truncated Levy random walks as shown in the paper. The 
nonlinear Langevin equation governing the particle motion was solved numerically using an order 1.5 
strong stochastic Runge-Kutta method and the obtained numerical data were employed to calculate the 
geometric mean of the particle displacement during a certain time interval and to construct its distribution 
function. It is demonstrated that the time dependence of the geometric mean comprises three fragments 
following one another as the time scale increases that can be categorized as the ballistic regime, the Levy 
type regime (superballistic, quasiballistic, or superdiffusive one), and the standard motion of Brownian 
particles. For the intermediate Levy type part the distribution of the particle displacement is found to be 
of the generalized Cauchy form with cutoff. Besides, the properties of the random walks at hand are shown 
to be determined mainly by a certain ratio of the friction coefficient and the noise intensity rather then 
their characteristics individually. 
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1 Introduction 

Levy random walks and Levy flights are met in a large 
variety of systems different in nature (see, e.g. [HHIEIUMSI 
6,7,8]). In the strict mathematical sense, such Markovian 
stochastic processes are characterized by divergence of the 
second moment of walker displacement x(t) for any time 
scale t, i.e., ([x(t)] 2 ) — > oo, which is caused by power-law 
asymptotics of the distribution function P(x,t). For ex- 
ample, in the ID case it is P(x,t) ~ Dt/x 1+a for x 3> 
(Di) 1 /", where D is some constant and the exponent < 
a < 2 specifies the superdiffusion law usually written 
in a symbolic form as ([x(t)] 2 ^ oc t" with the exponent 
/? = 2/ a > 1. Naturally, in the reality this divergence 
should be suppressed, for example, by the finite size of a 
system at hand. Thereby the given asymptotics can hold 
only within the frameworks of certain spatial and tem- 
poral scales bounded from below and above. The notion 
of truncated Levy random walks (flights) takes into ac- 
count the existence of these boundaries. To allow for such 
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cutoff effects explicitly several approaches have been pro- 
posed. These include continuous-time random walks gov- 
erned by the coupled spatial-temporal memory with finite 
moments [S?lllO|. coupled continuous-time random walks 
with bounded variations in the velocity fluctuations |11) . 
a direct cutoff in the Levy noise [T2j[T3] , Levy flights con- 
fined by external potentials [2], as well as Levy flights 
damped by dissipative nonlinearity [To] . 

Papers [TdIITTIITS] developed a new approach to de- 
scribing Levy random walks based on continuous Marko- 
vian nonlinear stochastic processes. For example, in the 
ID case it deals with random motion of a particle along 
the re-axis whose velocity v obeys the following stochastic 
differential equation 

= —\v + \j2r(v\ + v 2 ) * (la) 

written in the Hanggi-Klimontovich form, which is indi- 
cated with the multiplication symbol * in the product of 
the white noise £ (t) and the intensity yj2r{v 2 + v 2 ) of the 
Langevin forces depending itself on the particle velocity v. 
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Here r is a certain "microscopic" time scale of the parti- 
cle dynamics, A is a friction coefficient, and the parameter 
v a actually quantifies the intensity of "additive" noise. In- 
deed, a time patterns {v(t)} with the same statistics can 
be generated using also the following stochastic differen- 
tial equation 

dv 



T— = — 



dt 



on anomalous relaxation phenomena and introduced gen- 
eralized Cauchy processes appealing to a power-law decay 
in correlation functions. For example, the Cauchy distri- 
bution meets this interpretation when zt cx t@/ 2 : because 
in this case -Pc(O) oc t~P . Konno et al. [TW2"T] related a 
generalization of Cauchy processes to a power-law asymp- 
totics (maybe with some cutoff) of the distribution func- 
X + \f2r^ m {t) * v + V2rv a £ a (t) (lb) tion P{z) cx |z|~' 3 ' when the analyzed random variable z 
J takes large values, \z\ z t . For the Cauchy distribution 

the equality ft' = 2 holds. Below the term of the gener- 
alized Cauchy processes will be understood in the latter 
sense. 

Let us discuss a modification of the Langevin equa- 
tion (la I within the introduction of a nonlinear friction, 
i.e., a friction coefficient X(v) depending on the particle 
velocity v as well as the intensity G(v) of the Langevin 
random forces depending on the particle velocity v in a 
more complex way then it is accepted in model (la I. i.e., 
within the replacement 



explicitly containing additive, £ ffl (t), an d multiplicative, 
£ m (i), noise [TS]. It should be noted that equation (lb I 



has been much employed to study stochastic behavior of 
various nonequilibrium systems, in particular, lasers |20j . 
on-off intermittency [5T], economic activity [52], passive 
scalar field advected by fluid [23], etc. Besides, appealing 
to qualitative arguments and a numerical example Sak- 
aguchi [23] demonstrated that this equation can give rise 
to Levy flights in the space x on time scales t ^> r if the 
random variable v is regarded as the velocity v — dx/dt 
of a Brownian particle. 

on one hand, 



Random walks generated by model ( la 



can be treated as continuous trajectories of particle mo- 
tion and the corresponding probability density function 
P(x,v,t) obeys the standard Fokkcr-Planck equation. It 
admits a conventional generalization taking into account 
possible medium heterogeneities as well as the existence 
of system boundaries via the appropriate boundary condi- 
tions. On the other hand, on time scales f>r these ran- 
dom walks exhibit the characteristic properties, namely, 
the scaling law ([x(t)} 2 } cx t@ and the power-law asymp- 
totics of P(x, v, t) as x — > oo, that enable us to categorize 
them as Levy random walks. It has been proved analyt- 
ically for the superdiffusive regime (1 < < 2) [T7] and 
verified numerically for the quasiballistic and superballis- 
tic regimes (j3 > 2) [TH]. So the given approach seems to 
make it possible to attack the yet unsolved problem of the 
formulation of accurate boundary conditions for the frac- 
tional Fokker-Planck equations describing Levy processes 
in finite domains and heterogeneous media. 

The purpose of the present paper is to generalize the 
developed approach in such a way that explicitly allows 
for not only the small scale cutoff of Levy random walks 
determined actually by the parameter t in the model at 
hand but also a large scale cutoff. 



Generalized Cauchy processes 

We employ generalized Cauchy processes for constructing 
a model for truncated Levy random walks. There are sev- 
eral ways of introducing such stochastic processes starting 
from the Cauchy distribution 



Pc(z) cx 



1 



of a random variable z, where Zt is a certain parameter 
which, in principle, can depend on time t if, e.g., {z(t)} is 
some stochastic process. In particular, Lim and Li [35] fol- 
lowing Gneiting and Schlather [26 focused their attention 



A — > X(v) and 



G(v) 



(2) 



It should be pointed out that the introduction of nonlinear 
friction into a stochastic differential equation with addi- 
tive Levy noise gives rise to truncated Levy flights |15j as 
well as the Langevin equation (la) modified in the same 



way describes the generalized Cauchy process with a large 
scale cutoff [27] . 

Giving consideration to each of the two replacements 
has its own reason. First, the velocity dependence of fric- 
tion coefficient is widely used in describing nonlinear dis- 
sipation processes in nonequilibrium systems and for the 
systems where the intensity of dissipation grows with par- 
ticle velocity the Ansatz 



X(v) = X Q + X 2 v 2 



(3) 



with the coefficients Ao, A2 > is a simple and rather nat- 
ural model for these phenomena (see, e.g., [25]). Moreover, 
in describing motion of biological objects models assum- 



) . Appealing to 



ing X(v) cx \v\P are also met (see, e.g. 
Brownian motion governed by nonlinear friction and Levy 
noise |15| we may expect that a model similar to (la) with 
the friction coefficient ^ also generates truncated Levy 
random walks on large time scales. 

Second, the cutoff effect can be also due to the in- 
tensity G(v) of the Langevin random forces reaching its 
saturation as the velocity v increases. By way of example, 
let us mention movement of animals in searching for food 
resources, mates, den sites, etc. It has been found out that 
when animals have no information about the targets the 
resulting movement patterns of many species are of frac- 
tal, i.e., scale-free structure at least within multiple scales 
and can be described in terms of Levy flights or trun- 
cated once (see, e.g., [311] )■ The complexity of the move- 
ment patterns is due to combination of various "behav- 
ioral modes" that change over time and, thus, this animal 
movement can be regarded as a composition of Gaussian 
random walks with the intensity changing in time |30j (see 
also the following discussion [3"Tll3"2"] ). To quantify the in- 
tensity of the current behavioral mode the magnitude |v| 
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V a V c 

particle velocity, V 



some common "reference point" in comparing the cutoff 
effects caused by these mechanisms we will consider the 
models characterized by the same stationary distribution 
function. In other words, the parameters Ao, A2, and e will 
be set equal to such values that the ratio \{v)/G 2 (v) be 
the same function in both the cases. In addition, to single 
out the key points of the models at hand let us convert to 
the dimensionless variables, namely, 

t — > tr , x — > (v a r)x , v — > v a v . (8) 

Below all the results will be presented using these dimen- 
sionless time t, spatial coordinate x, and velocity v. 



Fig. 1. The characteristic dependence of the noise intensity 
on the particle velocity that mimics the animal behavior in 
foraging. It was used in constructing Ansatz Q. 



of the animal velocity v can be used as a natural param- 
eter that aggregates in itself the ultimate stimuli to the 
given behavior. So in mimicking the animal motion pat- 
terns based on stochastic differential equations similar to 
( la) the noise intensity G(v) may be approximated by the 
Ansatz 



G(v) = 



1 



ev* 



(4) 



It takes into account the low threshold in the animal per- 
ception of the optimal behavior as well as the animal 
bounded capacity of intensifying the search process quan- 
tified by the parameter e such that ev 2 <C 1. Figure [l] 
illustrates dependence Q. 

For the stochastic process governed by equation ( |Ta[ ) 
within replacements ^ the probability distribution func- 
tion P(v, t) obeys the following Fokker-Planck equation 



dP 

1h 



d_ 
dv 



G 2 (v) 



op 

dv 



v\(v)P 



(5) 



Its steady state solution, i.e., the stationary distribution 
P st (v) is of the form 



P st ( v ) = Cexp 



u\(u) 



du 



(0) 



where C is the normalization constant specified by the 
equality 



1 

C 



2 / cxp 





u\(u) 
CP(u) 



du 



dv 



(7) 



2 Model 

Random walks {x(t)} of a wandering particle in the ID 
space x S M are under consideration. The dynamics of 
its velocity v ~ dx/dt is assumed to be governed by the 
following Langevin equation of the Hanggi-Klimontovich 
type 

^ = -(a + l)vk(v) + V2g(v)*at)- (9) 

Here is the white Gaussian noise with the correlation 
function 

mm)=s(t-f), ao) 

the value < a < 2 is a model parameter and k(v) > 0, 
g{y) > are certain positive definite smooth functions of 
the argument v 2 such that 



fc(0) = 1 and 5(0) = 1 



(11) 



The numerical coefficient y/2 has been introduced into Q 
for the sake of convenience. In what follows the functions 
k(v) > and g(v) > will be referred to as the friction 
coefficient and the noise intensity, respectively, or just ki- 
netic coefficients with respect to both of them. 

When the kinetic coefficients are specified by the ex- 
pressions 



k (v) = 1 and g Q (v) = \f\ + 



(12) 



the stationary distribution function ([6]) which in the given 
case takes the form 



1 



(l+a) 



(13) 



stemming directly from the normalization of the distribu- 
tion function P st (v) to unity. 

In what follows we intend to analyze the two possible 
mechanisms of the cutoff individually. Namely, two alter- 
native cases will be studied employing separately the for- 
mer replacement of ([2| or the latter one. In order to have 



exhibits power-law asymptotics and its second moment di- 
verges for a < 2. Here B(-, •) is the Beta function. Under 
such conditions stochastic motion of this particle can be 
regarded as Levy random walks on time scales t ^> 1 [TBI 
H7|IT8]. Moreover, for a G (0, 1) the generated Levy ran- 
dom walks are of the superballistic type [TS] , whereas for 
a E (1, 2) they are of the superdiffusive type [T6lll7j . When 
a — 1 the stationary distribution function of the particle 
velocity is of the Cauchy form and the corresponding ran- 
dom walks can be treated as quasi-ballistic ones [18 . 
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As stated in Introduction, in the present paper the 
cutoff effects are related to deviation of the kinetic coef- 
ficients k(v) and g(v) from the ideal dependences ko(v) 
and go(v). To be specific we will make us of the following 
Ansatze 



k c (v) = 1 + 



and g c (v) — 



1 1 + v 2 
1 + v 2 /v 2 c 



(14) 



where the parameter v c 3> 1 is a certain critical value of 
the particle velocity when nonlinear effects responsible for 
the cutoff become essential. 

Below we will consider in detail two cases, 

(k(v) = k c (v) . TT fk(v) = k (v) 
case I < ) . ; ( and case II => I ) [ , , ; 
[9(v) = g {v) [g(v) = g c (v) 

which enables us, in particular, to elucidate the difference 
between the Levy type random works generated by the 
stochastic process Ety when either the friction coefficient 
or the noise intensity gives rise to the cutoff. The case 
k(v) — ko(v) and g(v) = go(v) will be also referred to as 
the unlimited Levy random walks. 

As intended, in the two cases the stationary distribu- 
tion function ([6| is of the same form 



Ce- 



ll 1 \ (!+») 



(15) 



(1 + v 2 ) 



with the normalization constant given by the expression 

, (16) 



C P = 



1 2- a (1 + a) (1 + a) 



2' 



2v 2 



2v 2 



where U(-,-,-) is the confluent hyp ergeometric function of 
the second kind. Distribution ( 15 ) can be regarded as a 



generalized Cauchy distribution with cutoff whose charac- 
teristics are studied in detail in Ref. |2"T1. 



3 Results of numerical simulation 

To analyze the statistical properties of random walks gen- 
erated by model Q it has been solved numerically us- 
ing the SRI2W1 algorithm of the Runge-Kutta methords 
for strong approximation of stochastic differential equa- 
tions of the ltd type with scalar noise [33]. This algo- 
rithm has the deterministic order 3.0 and the stochastic 
order 1.5. Mersenne Twister algorithm by Matsumoto and 
Nishimura [33] implemented in GNU Scientific Library 
1.14 }35j was used in generating random numbers. In in- 
tegration the time step was set equal to dt = 0.01, which 
gives stable results with respect to a decrease or increase 
of the time step by several times. The integration time 
was chosen to be equal to T = 10 8 in order to accumu- 
late enough statistics for analyzing the distribution tails 
based on individual trajectories, which enables us to avoid 
the problem of possible partial ergodicity of systems un- 
der consideration [18J . The created computer program was 



verified, in parti cula r, by reproducing the velocity distri- 
bution function ( |15[ ) numerically. 

The obtained numerical data were used in order to 
analyze two characteristics. The first one is the geometric 
mean x g (t) of the particle displacement 5tX during the 
time interval t: it is defined via the formula 



\n[x g (t)] = (\n(\S t x\)) 



(17) 



Here the angle brackets (...) denote the time averaging, 
i.e., partitioning of a generated trajectory of length T ^> t 
into fragments of duration t and then averaging over the 
obtained fragments. In the given model for the unlimited 
Levy random walks, i.e., for v c = oo the asymptotics of 
the geometric mean x g (t) as t — > oo can be written in the 
form [T7|[T5] 

x L g {t) = A g t 1 ' a , (18) 

where the coefficient 



r 



2- a 



r(a) 



exp 



(1-a) 



(19) 



r(. . .) is the gamma function, and 7 ~ 0.5772 is the Euler- 
Maschcroni constant. 

The second characteristics is the probability density 
P(p, t) of the particle displacement Stx during the time in- 
terval t measured in units of the corresponding geometric 
mean, i.e., the random variable p = S t x/x g (t). Again this 
distribution density was constructed based on partition- 
ing of a single sufficiently long trajectory. For the unlim- 
ited Levy random walks the asymptotics of the probability 
density is known, namely, it is of the form |18j 



Pl(p) 



2 /7TCV , 

sin ( — ) r(l + a)e 



7(a-l) 



1 



(20) 



as p — > 00 and does not depend on the time scale for t 1. 

Figure [2] shows the obtained results for the probabil- 
ity density P(p,t). Appealing to this figure we can state 
the following. When the cutoff region is rather distant 
from the origin in the velocity space, i.e., the parame- 
ter v c is rather large, model ^ descries random walks 
for which the particle displacement Stx obeys the general- 
ized Cauchy statistics with a certain cutoff. The solid line 
in the upper frame of Fig. [2] visualizes asymptotic ( 20 1 
corresponding to the unlimited Levy random walks. In 
this frame the dotted lines display the phenomenological 
Ansatz 



P(p,t)~P L (p) exp 



as p 



(21) 



where p c sa 3.2 x 10 2 and p c 3.3 x 10 for v c = 10 5 
and v c = 10 3 , respectively. The found magnitudes of the 
parameter p c validate the following construction of the 
relationship between the critical velocity v c and the char- 
acteristic scale A c determining the cutoff region Stx > A c 
in the space of the particle displacement {Stx}. 



Ihor Lubashevsky: Truncated Levy Random Walks and Generalized Cauchy Processes 



5 



10° 

10' 1 

o 

'■a io 2 -= 

c 

£ io- 4 

■o 

i 1 io- 5 

o 

* 10- 
10" 8 

10° 

I 10 " 1 

u 

I 10" 2 

J-J 

t/i 

10" 3 

I 10" 4 

n 
o 

a. 10 s 
10" 6 




case I: v c = 10000 
case II: v c = 10000 
case I: v c = 1000 
case II: v c = 1000 
case I: v c = 100 
case II: v c = 100 
asymptotics of LRW 
ansatz LRW x exp - 
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10 
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100 1000 



0.1 1 

particle displacement (in units of geometric mean) 
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♦ case I: v c = 1000 a = 1.0 

* case II: v c = 1000 t = 100 



case I: v c = 10 
case II: v c = 10 
Gaussian PDF 



case I: v c = 100 
case II: v c = 100 



5 10 15 20 

particle displacement (in units of geometric mean) 

Fig. 2. Probability density P(p, t) as a function of the particle 
displacement normalized to the the corresponding geometric 
mean, p{i) = 5tx/x g (t), for various values of the cutoff velocity 
v c . The used values of the parameter a and the duration t of 
the analyzed trajectory fragments are shown inside the frames. 



As it is known [16U17) , for the unlimited Levy random 
walks there is a relation between the particle displacement 
5tx attained within a given trajectory fragment of dura- 
tion f > 1 and the extreme fluctuations in the particle 
velocity within this fragment. Namely, the magnitude of 
5 t x is mainly determined by motion of the particle during 
the peak of the corresponding time pattern {v(t')} with 
the maximal amplitude i;J nax . Since the characteristic time 
scale of such extreme velocity fluctuations is about unity 



(in dimcnsionless units), the proportionality Stx 
with c ru 1 holds. Moreover, the estimate 



2 sin 



r 



(22) 



can be employed [161 and for a = 1 we have c = 2. In 
particular, the geometric mean x g (t) actually evaluates 
the characteristic amplitude of these extreme peaks gained 
by the particle velocity during the time interval t, i.e., 



civ. 



A*) mm- 



The cutoff effect in the statistics of the particle dis- 
placement has to become pronounced when the amplitude 
of extreme velocity fluctuations gets the cutoff region in 
the velocity space, ~ v c . It immediately enables us 

to write the desired relationship 



A c ~ cv c 



in the used notations it becomes 



P C (t) ■ Xg(t) - CV C 



(23a) 



(23b) 



According to the results to be presented below as well 
as obtained in Ref. [TH], the estimate x g (t) w 100 holds 
for the unlimited Levy random walks when the parameter 
a = 1 and the time scale t = 100. The found magnitudes 
for the quantity p c and the geometric mean x g together 
with the values of the para mete r v c used in the numeri- 
cal simulations fit estimate (23b) very well, justifying the 
given arguments. 

As the time scale t increases the intermediate region 
matching the Levy type asymptotics (20 1 should shrink 



and disappear completely when the equality x g (t) ~ v c if 
achieved. This conclusion is justified by the obtained nu- 
merical data illustrated in Fig. [2j It plots the probability 
density functions of the particle displacement that were 
constructed for the fix time scale t = 100 but different 
values of the parameter v c . Naturally, for x g (t) 3> v c the 
probability density function P(p) must be of the Gaussian 
form, which also is demonstrated in Fig.[2j the lower frag- 
ment. It should be pointed out that the crossover from the 
generalized Cauchy distribution with cutoff to the Gaus- 
sian distribution cannot be represented as a continuous 
transformation of Ansatz (21). In fact, the probability 
density P(p) for the intermediate values of the parame- 
ter v c — 100 « x g {t) for a = 1 and t — 100 looks like an 
exponential function. To make it clear the Ansatz 

P(p) = 1.2e-°- 74p 

is depicted by the dashed line in Fig.[2j the lower fragment. 

In addition, Figure [2] demonstrates the fact that statis- 
tics of the particle displacement depends rather weakly on 
the individual details of the kinetic coefficients k(v) and 
g(v), only their ratio k(v)/g 2 (v) matters to it. It could be 
explained if the statistics of the extreme velocity fluctua- 
tions is determined mainly by the stationary distribution 
function of the particle velocity, which is worthy of indi- 
vidual investigation. 

Figure [3] visualizes the time dependence of the geo- 
metric mean x g (t) for several values of the parameter a, 
namely, a = 0.5, 1.0, and 1.5. These values match the su- 
perballistic, quasiballistic, and superdiffusive type of the 
expected Levy scaling law x g (t) cx t x l a . The choice of the 
corresponding values of the cutoff velocity v c will be ex- 
plained below. 
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Fig. 3. The geometric mean of the particle displacements vs 
the duration of time steps in the partitioning of the random 
walks. The symbols □ and • depict the data points for cases I 
and II, respectively. 



Appealing to the data plotted in Fig. [3] we can single 
out three characteristic stages in the x g (i)-dependence. 
First, on small scales the time dependence x g (t) is found 
to be linear, as it must because of the strong velocity cor- 
relations for t < 1. Second, there is a certain intermedi- 
ate stage that can be classified as the region of the Levy 
scaling law. In fact, according to Ref. [TH] random walks 
governed by model Q with v c = oo exhibit the Levy type 
scaling law (18 1, i.e., x g (t) — Agt 1 / 01 , when the duration 
of the trajectory partitioning fragments exceeds a certain 



value about 100, i.e., t > 100 (in dimensionless units). 
So to make the intermediate Levy asymptotics feasible, 
in the numerical simulation the critical velocity v c was 
chosen to exceed the geometric mean x g (t)\t=wo for un- 
limited Levy random walks tenfold. The employed values 
v c meeting this requirement are noted in Fig. [3j In order 
to clarify the fact that in the cases under consideration 
the geometric mean x g (t) does exhibit the intermediate 
Levy asymptotics the function A g t x l a is also plotted in 
Fig. [3j The lines visualizing this function are seen to be 
practically tangents to the curves interpolating the numer- 
ical data. Finally, for large values of time scale the cutoff 
effect becomes substantial and, as it must, the region of 
the classical asymptotics x g (t) oc t 1 / 2 of Brownian motion 
arises. 

It should be pointed out that in both cases I and II 
actually the same geometric mean x g was obtained, which 
again justifies the previously drawn conclusion about the 
weak dependence of the given truncated Levy random 
walks on the individual details of the kinetic coefficients 
k(v) and g(v). 



4 Conclusion 

The paper is devote to the development of a continuous 
theory of Levy random walks that was recently initiated in 
Refs. pj5l[T71[T5] . Namely, a continuous Markovian model 
for truncated Levy random walks has been developed. It is 
based on the nonlinear stochastic differential equation Q 
governing the velocity v of a wondering particle. The cut- 
off effects in the distribution of the particle velocity and 
displacement during a certain time interval can be caused 
by the velocity dependence of the friction coefficient as 
well as the saturation in the growth of the Langevin force 
intensity with the particle velocity. Both the mechanisms 
are studied individually and to be able to compare them 
against each other the system parameters have been cho- 
sen in such a manner that the velocity distribution be 
the same in both the cases. The one-dimensional system 
was analyzed numerically employing a strong stochastic 
Runge-Kutta methor of the stochastic order 1.5 The 
obtained numerical data were used to calculate the geo- 
metric mean x g (t) of the particle displacement Stx during 
the time interval t and the probability density function 
P(p, t) of this displacement normalized to its geometric 
mean, p = 5tx/x g (t). 

It has been demonstrated that the two mechanisms do 
give rise to truncated Levy random walks on time scales 
exceeding substantially the "microscopic" time character- 
izing strong correlations in the velocity fluctuations. It is 
the case when the region v > v c , wherein the friction coef- 
ficient nonlinearity as well as the noise intensity saturation 
become substantial, is rather distant from the origin in the 
velocity space. As far as particular results are concerned, 
the following should be noted. 

First, the time dependence of the geometric mean x g (t) 
has been found to contain three characteristic asymptotics 
following one another as the time scale t grows. The initial 
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one matching the ballistic regime x g (t) cx t takes place for 
t < 1. The intermediate Levy scaling law x g (t) cx t x l a , 
which can represent the superballistic, quasiballistic, or 
superdiffusivc regimes depending on the parameter a G 
(0,2), is implemented for the time scales t > 100. On 
larger scales the cutoff effects become substantial and the 
x s (t)-dependence converts into the classical law of Brow- 
nian motion, x g (t) cx i 1 / 2 . 

Second, in the region of the Levy scaling law the con- 
structed probability density of the particle displacement 
5tx has been shown to be of the generalized Cauchy form 
with an exponential cutoff pronounced in the region S t x > 
A c . Using the numerical data the relationship A c ~ v c be- 
tween the cutoff scale in the space of particle displacement 
and the velocity v c characterizing the nonlincarity of the 
kinetic coefficients is justified. 

Finally, the properties of the given truncated Levy ran- 
dom walks have turned out to depend weakly on the in- 
dividual characteristics of the friction coefficient k(v) and 
the noise intensity g(v). Only the ratio k(v)/g 2 (v) spec- 
ifying directly the stationary distribution of the particle 
velocity matters to them. 
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